Quantum process reconstruction based on mutually unbiased basis 
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We study a quantum process reconstruction based on the use of mutually unbiased projectors 
(MUB-projectors) as input states for a D-dimensional quantum system, with D being a power of 
a prime number. This approach connects the results of quantum-state tomography using mutually 
unbiased bases (MUB) with the coefficients of a quantum process, expanded in terms of MUB- 
projectors. We also study the performance of the reconstruction scheme against random errors 
when measuring probabilities at the MUB-projectors. 
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I. INTRODUCTION 
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In physical implementations of quantum information pro- 
cessing, the ability for reconstructing the quantum process 
^^associated with a particular experiment is one of the most 
I l important tasks, which is mainly related to the determi- 
nation of the average fidelity of a quantum process in the 
Hilbert space of states. The full characterization of the im- 
perfection operations on a sequence of unitary transforma- 
tions in quantum devices, and the effects of decoherence on 
QQ the desired quantum evolution allow one to determine the 
fidelity of the quantum process. One of the more commonly 
used procedures is the so-called standard quantum process 
tomography (SQPT) [H-Qj where the dynamics of a quan- 
tum black box, described by a completely positive map £{p) 
T— I [IjIE], is reconstructed. In this method, for a D-dimensional 
^ quantum system, one is required to prepare an ensemble of 

•'H linearly independent input states {pk = |^fe)(V'fc|}fcLo 
subjecting each of them to the quantum process to be char- 
acterized, followed by standard quantum-state tomography 
■on the output states, £{pk)- The hnearity of £ relates the 
dynamical map to the experimental outputs via a linear sys- 
tem of equations which, after an inversion of a linear system 
of equations, allows to reconstruct the quantum-dynamical 
map. 

Several experiments have been implemented by consider- 
ing the above described reconstruction scheme, among them 
liquid-state NMR [ei], hght qubits 0,11, atoms in optical lat- 
tices tSl, cavity QED ^10,], trapped ions [ll| and solid-state 
qubits [12|. In all these experiments, the use of SQPT al- 
lows: the determination of the Kraus operators associated 
with the system's dynamics Il2l ; estimation of the presence 
of decoherence mechanisms [l3|; com put ing the experimen- 
tal fidelity of a quantum gate @, Isl. and, in this way. 



we can have an estimation of how close the real dynamics 
are to the theoretical prediction. In SQPT only single-body 
interactions are required for the reconstruction procedure. 

Recently, an alternative scheme for quantum process to- 
mography (QPT) which uses an ancillary system has been 
implemented, which is the so-called assisted process tomog- 
raphy (AAPT) [il-[l3. In this scheme, an isomorphism be- 
tween a quantum process and a quantum state in an extended 
Hilbert system is established jlTj]. An input state, defined 
on the joint Hilbert space, passes through the map £ and 
a projective measurement on the the output state is per- 
formed, thus reconstructing the quantum-dynamical map. 
These measurements can be implemented by using factor- 
ized measurement or joint measurements using mutually un- 
biased bases (MUBs). If a quantum map acts on a set of 
qubits, one ancilla is added to each qubit and joint mea- 
surements are performed on one system's qubit and on its 
ancilla. For this reason, for an iV-qubits system A^-ancillas 
are needed, and 5^ measurements must be carried out for 
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the reconstruction of the process 17j . Another scheme that 
makes use of MUBs is the selective and efficient estimation of 
the parameters that characterize a quantum process, using 
MUBs for averaging the fidelities of appropriately modified 
channels fH, [l^ . This scheme also needs an ancillary system 
for the determination of the off-diagonal terms of the process 
in its matrix representation. 

In this article, we provide a new scheme for QPT, which 
is based on mutually unbiased basis quantum-state tomog- 
raphy (MUB-tomography) over the output states [l^l- This 
procedure is valid for a Hilbert space with prime-power di- 
mension, i.e., D = with p and r being a prime number and 
an integer, respectively [2l[. For instance, it can be applied 
to a map acting on a single high dimensional system or on 
a multi-qubit system. In this scheme, two-body interactions 
are needed and do not require ancillary systems. 

This paper has been organized as follows: In section |TT] 
we present the reconstruction of the QPT using MUBs, or 
simply MUB-QPT, showing the expansion of the quantum- 



2 



dynamical map £ in terms of MUB-projectors. Furthermore, 
we derive a relation between the expansion of the quantum 
map and the measurement of MUB-tomography, via a lin- 
ear system of equations. In section Hill we present numerical 
simulations on the performance of the protocol against the 
presence of random errors in MUB-tomography. In the nu- 
merical simulation, we have considered quantum-dynamical 
maps describing a local decoherence channel and non-local 
quantum map. The article ends with a summary in Sec. IIVI 



II. QPT BASED ON MUB-TOMOGRAPHY 

For this purpose, let us consider the set of mutually unbi- 
ased projectors (MUB-projectors) for a I?-dimensional sys- 
tem, whose dimension is a prime or a prime power number 
{D = with p a prime number and r a positive integer), 

where a set of mutually unbiased bases \ipm^) exist [2l] - [26| . 
For a ZJ-dimensional system it has been found that the max- 
imum number of MUBs cannot be g reater than D -I- 1 and 
this limit is reached if D is prime [22l | or power of prime (2]| . 
Recently, it has been suggested that the existence of MUBs 
for other dimensions may well be related to the non-existence 
of finite projective planes of certain orders [27[ or with the 
problem of mutually orthogonal Latin squares in combina- 
torics [28i]. 

The MUB-projectors are given by: 



nl^^ = I^L^^)(^L^^I, m=l,...,A 7 = 0,..,A (1) 

where 7 labels one of the D + 1 families of MUBs and m 
denotes one of the D orthogonal states in this family. These 
projectors satisfy the following relation: 

Tr [v^7?Vlf^) = hiSrnn + ^ ' ^ m) ■ (2) 

Besides, these operators define a complete set of projection 
measurements, i.e., X]m=i^™^ = 1 (1 denotes the identity). 
The measured set of probabilities p^m — com- 
pletely determines the density operator of the system, so 
that 



D D 
m— 1 7—0 



(3) 



In a previous work [20j, we have shown an optimal protocol 
to obtain the probabilities p-ym, for the case of having a multi- 
qubit system. The reconstruction scheme is based on mini- 
mizing the number of conditional logic gates (two-particle 
quantum operations) used in quantum-state tomography, 
due to the fact that all the MUBs contain non-factorizable 
bases. In this work, our main goal is to perform the quan- 
tum process reconstruction in a fully symmetrical form using 
MUB for input states, output state reconstruction and to ex- 
pand the quantum dynamical map. The main motivation for 



doing this is the recent advancement in experimental imple- 
mentation of MUBs quantum state tomography. It has been 
experimentally demonstrated quantum state tomography of 
two-qubit polarization using MUBs [29]. They demonstrate 
an improved state estimation when comparing with standard 
reconstruction. Besides, an experimental implementation of 
quantum tomography of 7- and 8-dimensional quantum sys- 
tems has been reported (30| , where higher dimensional quan- 
tum systems are encoded using the propagation modes of sin- 
gle photons, and MUBs projections have been implemented 
using programmable spatial light modulators. 

To approach this problem, we begin by considering the 
general evolution of a _D-dimensional quantum system de- 
scribed by a completely-positive linear map £ (p) , expressed 
in the so-called operator-sum representation^ 



(4) 



where Ai are the Kraus operators of the system and satisfy 

To expand the Kraus operators of the system, we consider 
the overcomplete basis of {D'^ +D) MUB-projectors, Vm'' to 
expand the operators Ai, 



D D 



(5) 



Hence, the completely-positive linear map £ [p) can be ex- 
pressed in the following manner: 



a,/3— m,n— 1 



(6) 



where Xm£'' 



is the process matrix. If we de- 



mrf\ we can reconstruct 



termine the {D^ + D)^ quantities x 
the quantum process in terms of the set of MUB-projectors 
Vm"^ . Unless, the determination of the unital condition [sij, 

^ 1 and the features of MUB-tomography 
mainly related with the expansion of the density operator 
in terms of MUB-projectors, as is shown in Eq (jS]), implies 
that we need to use the overcomplete basis of inputs to ex- 
pand the Kraus operators and the density operator of the 
output states. In this way we can present this scheme as 
an extension of the results obtained previously for quantum- 
state tomography using MUB. On the other hand, due to the 
complete set condition for each family of MUB (labelled by 

a), the action of the quantum map over any input state Vi'^\ 
according to Eq. ^ and the linearity of the super-operator. 



E 



AAl - E ^ (^™^ 



(7) 



Then, the action of the map on the remaining element T'l'*'^ 
(m ^ /) could be calculated by using the results of the super- 
operator £ acting over the D — 1 elements, J2m=ii ^ ^) ■ 
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But the quantum dynamical map is unknown, and we can't 
predict the value of J2i ^i^l- other words, we don't know 
if the map is unital or not (sTj . For this reason, we need 

to determine the output state £ ^'P/^^^ by quantum-state 

tomography. To summarize, for the reconstruction of the 
quantum process we need to prepare {D^ + D) input states 

Vl'^'^ that correspond to the set of MUB-projectors. These 
states are sent through the quantum black-box and, after 
that, we carry out quantum-state tomography on the output 
states £{vl'^'^), based on mutually unbiased measurements 
[20]. We emphasize that the set of {D^ + D) MUB-projectors 
used to expand the Kraus operators is the same set used as 
input states for the quantum process. 

We can expand any density operator using the overcom- 
plete set of MUBs [20| . For this reason, the dynamical map 

associated with each input state, £{vl'^^) can be expressed 
in the following way according to Eq. Q: 



e=0 5=1 



(8) 



are linearly independent. This is not a problem, because the 
inversion of the linear system is guaranteed in terms of the 
generalized inverse, as we show in Sec. (jllip . even when only 
D"^ Eqs. are linearly independent. Also, there is the pos- 
sibility to use only MUB-projectors and the problem of 
QPT is solved, but in that case the extension of our previous 
result, shown in Ref. [13] related to the reconstruction of 
quantum process, is not clear and that is the reason to make 
a protocol based on an overcomplete basis set. 

Now, for simplicity we assume that the probabilities 
p]^/^ and the Xmn^ coefficients can be expressed as vec- 
tors {p and X, respectively) and that the values given by 

Tr can be arranged in a (D^ -f Df~ 

dimensional complex square matrix /3. 

Let us consider the inversion of equation ([TT|) . If we define 



Tr 



(7,0 

Pr)s , 



Pap\lp it follows that, for each 



_^(a,/3) nmn,ls 



(12) 



where pjq' = Tr{'Pq^'^ £{'P[^'')) are the probabilities associated 

with each Vq^^ projector after the dynamical map acts on 
the input state. These probabilities must be experimentally 
determined by using MUB-tomography. 

Also, we can write the same dynamical map using Eq. 

^i-Pl'')- E E xi^^n'^-Pt^VlM'^. (9) 
a,p=0 m,n=l 

Combining equations (|5]) and ^ we obtain 

j:f:pi,'H''-^- E E xi^n^'vL-'viM'^. (10) 

c—Oq—l a.P—Orn,n—l 

We make use of property 1^, i.e., we apply the vi^^ projector 
on the above expression and take the trace on the whole 
expression obtaining 

P'r//'=i: E xL";f (T'i-M^M^Pi'')) . (11) 

a,P=0 m,n=l 

We can observe from the equation ([5]) that for each input 
state vl'*'' we can obtain a set of £'(13-1-1) probabilities from 
the state tomography. Then, with the experimental results 
of all the input states, we get a set of {D^ + D)^ probabili- 
ties. Using this set, jJ^s''"*, we write a linear system of equa- 
tions pT|) with all the terms of the Xmn'^ process matrix and 
the coefficients associated with the corresponding probabil- 
ity of the set, given by Tr (^'l^^T'/^^'Pi^Vi''^) . Clearly, as 
is shown in Eq. ([7]), not all the Eqs. that arise from (fTTI) 



namely, p = /5x.To make a proper inversion of Eq. 
us consider k, the generalized inverse of matrix P (s! 
fulfills: 



nmn,ls nrt,xy ^v.ujO 

r^t.x.y 
fi v,uj ,9 



,mn,ls 



ni, let 
, which 

(13) 



Hence, Xmn^ is defined as, 



Amn 



E 

•y.rj.l^s 



ran.lsi^rjs 



(14) 



In this manner, the problem of QPT can be solved by finding 
the generalized inverse k of the complex matrix j3, and finally 
solving X = i^P- 

In general, the determination of probabilities in tomo- 
graphic reconstruction of output states is subject of exper- 



imental errors (pijl'''* — t- Pi;, 
denotes the error for each output state), which can give ori- 
gin to unphysical process matrix Xmrf^ when inverting the 
linear system, Eq. (IT2l) . This problem can be overcome in 
a similar way than in the quantum-state tomography Q by 
minimizing a deviation function f{t), where i is a general 
parametrization of the physical matrix process coefficients 
(t) [1, [l3- In our case, this function can be written as: 



(7:0 

P-qs 



where e 



(7,0 

r)s 



fit) = 



D D 



/ J / J /^mn \ } Amn 



a,/3— m,n 



(15) 



11\ 

1,1 



it) 



D D 



E ExL"^H*>K^"''^-'^' 



(7)^ 



Q,/3=0 rn.n 
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TABLE I: Five sets of three operators defining the MUB set for 
L> = 4 



OL 


Operators 





a -A 






1 








2 


ayi 


lay 


GyOy 


3 


d'xd'y 




CfyBz 


4 


ayCTx 


CfzCfy 





The second term in the right hand side is obtained by sum- 
ming on index s in Eq. (|12p and x'f^ are the Lagrange mul- 
tiphers corresponding to each input state. The parametriza- 
tion t is chosen in a such a way of preserving both the her- 
miticity and the positivity of the map. 

The advantage of mutually unbiased basis quantum-state 
tomography (MUB-Tomography) is related with the ability 
of both to maximize information extraction per measurement 
and to minimize redundancy. Hence, the estimation of the 
output state has several advantages with the use of MUB- 
Tomography. For instance, there is no need to invert a linear 
system of equations, as in standard quantum-state tomogra- 
phy j29i . iSOi] . Our main purpose is to extend this advantage 
to quantum-process reconstruction by having a fully sym- 
metrical form using MUBs, namely, using MUBs for input 
states, for output state reconstruction and for expanding the 
quantum dynamical map. 

In the next section, we give an explicit example of this 
protocol for local and non-local processes in Z? = 4 and show 
how the performance of the protocol is affected when random 
errors are introduced into the obtained probabilities from 
quantum-state tomography, for local and non-local channels. 



III. EXAMPLE IN D = 4 



which arise from the presence of noise in the estimation of 
the probabilities p^/\ in each quantum process considered. 

The Kraus operators that represent local decoherence op- 
erations we denoted as Aij. They act individually at each 
qubit (labelled by (i, j)), so Aij — Ai ® Aj and satisfy 
AijAlj = i. These local processes, for a single qubit 
operation, have the following representation: 



Depolarizing Channel 


Amplitude Damping 


A2 = y/J^x 
A3 = y/l^v 
^4 - yf 





The input states evolve in these local channels as £{Vi ) = 

AijVi'^'' {Aij)\ We have chosen these channels because 

none of the input states Vm'' are an eigenstate of the Aij 

operators and, for this reason, all the states will be affected 

by these decoherence mechanisms. This assumption will be 

not accomplished if, for example, we choose the bit-phase 

flip channel [if. In this case, the operators Aij associated 

with the channel are those in the third row of table HI This 

(2) 

implies that the input states Vm will not be affected by this 
local decoherence mechanism, because they are eigenstates 
of all the operators that belong to the set a = 2 and, in this 
way, the set of three projectors Vm'' will be a decoherence- 
free subspace (DFS) for the bit-phase flip channel. We avoid 
this by choosing the amplitude damping and the depolarizing 
channel, because none of the 20 MUB-projectors obtained 
from the table |T] constitute a DFS. 

On the other hand, the considered non-local operation is 
the controlled-NOT gate, that can be written as a coherent 
sum: 

UcNOT = \{i®i + i®ax+(J^®i-d:,®ax) ■ (16) 



Here we apply the above described scheme for the pro- 
cess reconstruction in the case of a four-dimensional Hilbcrt 
space, I? = 4, which could be associated to a two-qubit sys- 
tem. There are several approaches for the construction of 
mutually unbiased bases for this bipartite quantum system 
[2l| - |26j , which satisfies the prime-power condition for its con- 
struction, 2^ = 4. In this case, we will use the two-qubit 
operator set of table HI such that the common eigenstates 
constitute a MUB. Then, we can obtain 20 MUB-projectors 
from Table U (four from each row). 

The set of projectors v'f^ (7 = 0, .., 4 and Z = 1, .., 4) will 
be the input states for testing the accuracy of this procedure 
for three fixed quantum processes; amplitude damping and 
depolarizing channel, as local decoherence mechanisms [ij, 
and a controlled-NOT gate, Ucnot- Our purpose is to show 
the effects on the process matrix x associated with each map. 



This process evolves the input states according to £{Vi ) = 

{Ucnot)'Pi'^^ (Ucnot)^ ■ We apply this set of operations 

to the whole set of input states vl''\ as is shown in Eq. 

(H]). Then, with the output states obtained, £{vl'^^), we 
calculate the corresponding probabilities associated with this 

dynamical map, p^V^ = Tr{£{v'l''^)V'f'^). 

To have an estimation of the quality of this procedure, we 
have considered an error model for the probabilities p^/^ , be- 
cause in the experiments we always have to deal with small 
errors associated with these values. Then, our noisy proba- 
bility p^l''^ is obtained by adding to the calculated probabil- 
ity a random value Q (uniformly distributed between and 
1) times a constant error parameter /i: 

4"/^=H]''^+K^I''\ (17) 
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This quantity is defined as, 



F 



\ Xmn 



rr((xL"/^)2) 



(19) 



0.15 



FIG. 1: (Color online) Process fidelity F, for local and non-local quan- 
tum operations, when the process matrix x is determined via noisy 
measurements in the MUB-Tomography, with error parameter fi. Here, 
the fidelity in the determination of the process matrix associated with 
the non-local operation Ucnot is plotted with blue dashed-dot line. 
For the local channels, amplitude damping with 7 = 0.4 is plotted with 
red dotted line and the depolarizing channel with p = 0.1 corresponds 
to the black dashed line. 



where e^T''* 



(7,0 



This allows us to simulate a noisy 



measurement of the output states £{V^[''^), for the three dif- 
ferent quantum operations considered. For each output state 
we normalize the new probabilities obtained accord- 
ing to Eq. (|17p . for each measurement basis of the MUB- 



1. With 



Tomography (labelled by 77), that is, X)sP') 
these values, we can calculate a noisy process matrix asso- 
ciated with the quantum channels considered, Xrn,n \ which 
is obtained in an experiment with non-ideal measurements 
of the output states. We obtain this new process matrix by 
using Eq. (fT4|) with the noisy probabilities on the right side. 



(18) 



For this purpose, we perform a numerical simulation with 
100 steps for every value of /i, each one representing the error 
model for the probabilities defined by Eq. (fT7)) . i.e., it adds a 
set of random values into the probabilities of detection. The 
dimensionlcss error parameter ^ takes values between 0.01 
and 0.15, with increments of 0.01. We will then obtain then 
a set of 100 noisy process matrices Xmn\ for each value of 



The process fidelity F [s!, 'ss'l is a good way to estimate the 
quality in the determination of the dynamics with this basic 



error model, comparing the ideal process matrix Xmrf^ with 

Ymn^'' obtained in the simulation. 



each one of noisy matrices xJi 



and take values between and 1. Figure [T] shows how the 
mean fidelity on the determination of the process matrix de- 
creases when the error in the probabilities grows, as we ex- 
pected. The mean fidelity is calculated using all the noisy 
matrices calculated Xmrf "^ obtained in the simulation. Be- 
cause of the introduction of the random errors, the precision 
in the estimation of the gate decays, but in the case of non- 
local gates, this behavior is more notorious. To explain this 
issue, we take one entangled input state and then apply the 
noisy non-local operation. Then, we calculate how the con- 
currence of the state changes when we introduce the random 
error. We observe that the concurrence decays in the same 
way as the process fidelity for that operation. Hence, the 
quantum correlations that the non-local operation induces 
are destroyed by the introduction of the random error, mak- 
ing the process fidelity F decay more strongly than do the 
local ones. 

The measurements in any set of MUBs contain non- 
factorizable bases, we require application of nonlocal gate 
operations. These operations currently can not be performed 
with unit fidelity. In our previous work (20| . we defined the 
physical complexity, for each set of bases, as the number of 
CNOT gates needed for implementing projective measure- 
ments. Hence, for quantifying the amount of resources for 
implementing the QPT scheme, we use the physical complex- 
ity of the chosen set of MUBs. 

In the case of MUB-Tomography, for a _D-dimensional 
quantum system, the physical complexity for each set of 
MUB is denoted as Ca {a = 0,...,D). Then, for imple- 
menting projection measurements in (D + 1) bases, we need 
C = J2a=o CNOT Gates. Therefore, for each input state 
T'l'*'^ we need to perform C non-local operations to carry out 
MUB-Tomography over the output state £ (vl'^^^ ■ But, to 

prepare each input state (MUB-projectors) we also require 
a number of CNOT gates, given by the complexity of each 
MUB Ca- As D input states are prepared from each MUB; 
therefore the total number of CNOT gates required will be 
{DCa)- Summing over the {D + 1) MUBs, then we need 
{DC) non-local operations to prepare all the input states. 
Finally, we require a total number of (DC^) CNOT gates for 
the whole QPT protocol. 

Here, we assumed that the local operations are performed 
with unit fidelity. For a given dimension D of the Hilbert 
space, only certain factorization structures are admitted [23| , 
each one of them having a different complexity. Then, we 
need to minimize C in order to carry out a lesser number of 
CNOT Gates. This is discussed in more detail in [201 ■ 
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IV. SUMMARY 

The choice of the set of MUB-projectors both as input 
states and as an operational basis, and the calculation of 
pi^/^ by the experimental data of MUB-Tomography, per- 
mit us to reconstruct an unknown quantum process, deter- 
mining the process matrix Xmrf\ improving the way that 
standard QPT relates to the experimental outcomes with the 
parameters of the process matrix. We also can recover the 
Kraus operators for the system A^, by the diagonalization 
of the process matrix. The simulations suggest good process 
fidelity, even when the values of the random errors /^C^I''' 

are similar to the theoretically calculated probabilities p^^/\ 
according to Eq. p7)) . making this protocol robust against 
the presence of noise when determining probabilities. This 
protocol is valid in any Hilbert space with prime power di- 
mension and where two body interactions are available, due 
to the fact that MUB-Tomography for all input states re- 



quires (DC^) CNOT gates for its physical implementation, 
where C is the physical complexity of the chosen MUBs. 
The case of not completely-positive maps [33|, the perfor- 
mance of this protocol compared with other procedures of 
QPT discussed in Sec. U in terms of the inaccuracy of esti- 
mated probabilities in the experiments, and the dimension of 
the considered Hilbert space, and how these features involve 
the estimation of the quantum process, are now left as o pen 
questions. A similar experimental setup used in Ref. [30| 
can be used for a proof of principle implementation of this 
reconstruction scheme, where an extra spatial light modula- 
tor along the propagation path of single photons can be used 
for implementing different processes in a controlled way. 
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